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ABSTRACT 

A comparison of three methods to deduce a property of element potentials that allows to write the 
Gibbs energy expression for systems at chemical equilibrium is presented. These methods are: (1) 
the algebraic procedure widely used in reactive distillation modelling to obtain the transformed 
composition variables, (2) the classical non-stoichiometric formulation of chemical equilibrium 
problems, and (3) the complete Legendre transformation of the Gibbs energy with respect to 
element composition variables. These three techniques lead to the same result in most cases, i.e. 
when elements are chosen among the system components. Because of its consistency with the rest 
of the Thermodynamics formalization, the latter gives deeper insight into the transformed 
composition variables and the auxiliary function of the Lagrange multipliers method. The 
element composition variables are the conjugate variables of the element potentials. The auxiliary 
function of the Lagrange multipliers method is the complete Legendre transform of the Gibbs 
energy with respect to element composition variables, and its value is zero for equilibrium states. 
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1. INTRODUCTION 

In reactive mixtures, changes in composition variables are not independent. 
Thus, for reactive closed systems, the differential form of the fundamental 
equation takes a simpler form than the one for open non-reactive systems 
(Luetich 2002). It has been demonstrated that this simplification is due to the 
fact that such variables can be transformed in others which do not change when 
the equilibrium is displaced modifying the state variables. 

An algebraic method to carry out this transformation was introduced by Barbosa 
& Doherty (1987) with the purpose of transferring to reactive distillation 
techniques of non-reactive distillation modelling. 

Perez Cisneros et al. (1997) treated the same systems by means of the classical 
non-stoichiometric formulation (Smith & Missen 1991), showing that this 
method makes it possible to use composition variables which do not advance even 
in the case of reactions such as esterification, i.e. when not all the elements exist 
free in the mixture. 

At the same time, Alberty (1997) presented a method based on the use of 
Legendre transformations to treat systems containing a component with 
constant chemical potential. His approach was based on previous works related 
to the reduction of the number of composition variables in systems containing 
pseudoisomers (Alberty & Oppenheim 1990). 

In this work, it is demonstrated that using complete Legendre transformations of 
the Gibbs energy with respect to the element composition variables allows for a 
deeper interpretation of the physical sense of the transformed composition 
variables (Luetich 1999) and the auxiliary function of the Lagrange multipliers 
method (Luetich 2001). 

A comparison of the three techniques is done by applying them to the solution of 
a problem for which all lead to the same result. 

This paper is the second of four related works written to clarify the concept of 
microscopic reversibility. With this purpose, the emphasis here is laid on the 
information contained in the different forms of the fundamental equation for 
reactive systems. 
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2. WRITING THE GIBBS ENERGY EXPRESSION 

OUTLINE OF THE PROBLEM 

Given the reactive system studied by Brandani & DiGiacomo (1984) 

I formaldehyde + water ^ methylene glycol 
( methylene glycol ^ water + polyoxymethylene 

the simplest expression for the Gibbs energy is to be found. 
SOLUTIONS 

■ LINEAR COMBINATIONS 

To simplify the notation, generic components A, B, C, and D, instead of 
formaldehyde, water, methylene glycol and polyoxymethylene, respectively, will 
be used. 

, A+B^C 

{ 2 C^B + D <2 ' 2) 

The Gibbs energy of system (2. 2), 

m 

G = ( jU A /Ub He Mo)(n A n B n c n D ) , (2. 3) 

at equilibrium, is subject to the restrictions 
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(//a Hb He Hn) , „ =(0 0). (2.4) 



(This property accounts for the use of the adjective "chemical" for the magnitude 

H-) 

Then, only two chemical potentials are independent. 
(1 \ /-l 2 ^ 

(l o u-i 2 y 1 

(He Hb) = -(Ha ^)[ 1 _i II o -l) (2 ' 6) 

C and D are the so-called reference components (Barbosa & Doherty 1987). 
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On the other hand, associating terms in the matricial expression of the Gibbs 
energy, 

G = ( yU A fi B ) ( n A n B ) T + ( yu c /u B )(n c n B ) T . (2. 7) 

The following step is the most important one since, by means of it, two variables 
are eliminated: 

„ ( Ha + nc + 2 n B \ 

G = (ha fi B )\ A \. (2.8) 

This relation allows to define the transformed composition variables (h) . 

(Ha VB)( nA ) = G (2.9) 

\n B J 

In appendix A, it is shown that this procedure is equivalent to eliminate the 
degrees of advance of both reactions. 

■ LAGRANGE MULTIPLIERS 

The system can be viewed as constituted by "element" A and "element" B 
(physical approach): 

| A + B ^ AB 

l2AB^B + A 2 B (2 ' 10) 

Therefore, it is intended to minimize the function 

G = G{T,P,{n A n B n c n D )), (2.11) 

subject to restrictions 



(2. 12) 



which define the composition variables n* A and n* B . These variables arise by 
inspecting the number of "atoms" that conforms each molecule. 

The problem is just the classical minimisation one, with bound variables, solved 
by Lagrange using multipliers. 
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Let 



O = G- 



(A A A B ) 
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111 



n B 
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V ^D J 



(2. 13) 



be the corresponding auxiliary function [the Lagrangian, 0(n, A)]. 

Its derivatives with respect to mole numbers are 
( <90 ^\ 
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(2. 14) 



The matrix product should be transposed because 
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(2. 15) 



These derivatives become null -giving G its minimum value- when 
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(2. 16) 



Therefore, 



min G = (n A + nc + 2rij) n B + nc + nu ) 



(n* A n* B ) | A 1 = min G. 



A B J' 



(2. 17) 

(2. 18) 



The relation between both preceding methods is shown in appendix B. 



Trans. Luv. Acad. 2002, 1, 23-36 



28 



J. J. Luetich 



■ LEGENDRE TRANSFORMS 

For this reactive system, the Gibbs energy expression is 



G = (n A n B n c n B ) 



He 



(2. 19) 



At constant T and P, the equation which determines the change of this 
magnitude with composition variables is 



(dG) T p = ( yU A jU B jU C Ht>) d 



fn A \ 
n B 
n c 

V ^D J 



(2. 20) 



Then, the problem consists of finding out the minimum number of times in which 
the Gibbs energy should be transformed when variables are bound by chemical 
reactions stoichiometry so as to obtain its complete Legendre transform. But this 
number is just the number of elements. In other words, variables that remain 
constant during a reaction are the element composition variables. Thus, the 
transformation (2. 12) shoud be used. 



For this case, the complete Legendre transform of G is 



G-{n\ n* B ) 



^ n* 

= Lr t 



AB- 



(2.21) 



[A new symbol (£) is introduced here because it is not known yet which variables 
should go with the transformed composition variables. Variables £ are the 
element potentials.] 



Operating on equations (2. 21) and (2. 12) 
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(2. 22) 



(2. 23) 
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Since this equation holds for non null variations, we obtain the expressions for 
element and compound potentials: 



( ji A n B Mc Mb) = ( U £b £a + £b 2 £ A + & ). 



(2. 24) 



Therefore, 



G = (n* A n* B ) 



Ma 

Mb 



(2. 25) 



This is the same result obtained by means of the Lagrange multipliers method. 
Note that this technique -using integral equations- is the inverse of the linear 
combinations method (appendix C). 

From here on, to complete the deduction, it will be convenient to use the 
following notation: 

(Ma~£a Mb-£b /*c - & - £b Mb - 2 £ A - & ) = ( M I M\ M*c M*o )» ( 2 - 26) 
where fi* is the transformed potential. 

In addition, transformation (2. 21) defines a thermodynamic potential whose 
expression in terms of the original and the transformed composition variables, is 



G*^ = ( fie- Ha- Mb Mb -2 Ma- Mb) 



fn A \ 
n B 
n c 

V n B ) 



G* m = (0 m*c Mb) 



(n A \ 
n B 
n c 

V n B J 



(2. 27) 



(2. 28) 



which is zero at the equilibrium state. 



Thus, the differential form of this function is 
(dG^Tp = -(n\ n* B )i A j. 



(2. 29) 



Formulae, which relate the transformed and the original thermodynamic 
potentials and their derivatives with the transformed composition variables for 
equilibrium states, are shown in table 1. 
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Table 1 

f A+B^C 
Equilibrium conditions for the system < , in terms of the Gibbs and the transformed 

v 2 C ^ d + D 

Gibbs functions 



G = min G = n* A // A + n* B // B # 
(cIG)tp = /i A <^ A + y^B c^* B = 



G*^ = min G*^ = n* c n c + /u* B n B = 
{dG* m ) TP = -n\dn K -n\dn B = 



Both functions (G and G*^ ) have null differentials at equilibrium states, but only 
the transformed potential is zero. In the first case, the differential nullity is due 
to the invariability throughout the reaction of the transformed composition 
variables. In the second case, it is so, because the transformed composition 
variables are bound to potentials by a relationship analogous to the so-called 
Gibbs-Duhem equation. 

■ Generalization 

The general expression of the Gibbs energy for systems undergoing chemical 
reactions is 



G = G (m, n) = n T ft, 



(2. 30) 



where m is the vector of the state variables and n is the vector of the composition 
variables. 



m = (T, P) 



(2.31) 



Stoichiometric constraints reduce the number of composition variables: 

Cn = n*. (2.32) 

C is the conservation matrix. 

The intensive variables which go with the transformed composition variables in 
the shorten expression of the Gibbs energy can be related to chemical potentials. 



G = (n*) T <r 

t i* = f i-C T £= 



(2. 33) 

(2. 34) 



The general integral and differential forms of the transformed Gibbs function are 
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G* = {n *) T n; (2. 35) 

(dG*) m = -(n*) T c^. (2.36) 

From equations (2.34), the components of C, can be obtained. When replacing 
them in equation (2.36), they provide the expression of dG* according to the 
original chemical potentials (fi). The transformed potential results from 
equations (2.34). If they are replaced in equation (2.35), the integral expression 
of G* can be obtained. 

Otherwise, adding to equation (2. 23) the Gibbs-Duhem equation, 

(H*) T dn + iidfi = 0; (2.37) 

[d (G - (n*) T £)] m = (dG*) m = - (n*) T d£. (2. 38) 

Therefore, we obtain the equation: 

(dG) m = fdii*, (2.39) 

whose importance is apparent when studying chemical equilibrium problems 
within the frame of the physical approach (Luetich 2002). 

The transformed composition variables are the conjugate variables of the 
element potentials (Luetich 2001). 

Table 2 presents a summary of the three methods used to solve this problem. 
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Table 2 

One problem, three focuses 



THE PROBLEM 


ALTERNATIVE FORMS 


Given the reactive system 


Linear combinations 


< A+B^C 
120B + D 


Transform the mole numbers of C and D 


to put them as combinations of 


the simplest 


the ones of A and B in the expression of G. 


expression for the Gibbs 


G = n A ^i A + n B ^ B + n c /ic + ^d j«d 


energy is to be found. 


f fi A + /J. B ~ yU C = 

I -/i B + 2/i c -/id = 






Lagrange multipliers 




Using the element 




composition variables n* A and n* B , 




give an expression to the minimum 




of the function 




G = G(T,P,(n A n B n c n B )). 




Legendre transforms 




Find the conjugates of element composition 




variables to write the differential form of 




G = G(T,P,(n* A n* B )). 



3. CONCLUSIONS 

In the method introduced by Barbosa and Doherty (1987), the composition 
variables come after defining the way in which chemical potentials will be 
combined (equilibrium condition). On the contrary, in the Lagrange multipliers 
and the Legendre transforms methods, the selection of variables to be used in 
describing the system composition is made before obtaining the expressions for 
potentials. That is what allows us to find a simple expression for the 
thermodynamic function that defines the equilibrium states in the case of 
reactions such as esterification (Perez Cisneros et al. 1997), i.e when elements 
are always combined. 

When "elements" are chosen in such a way that information on the system 
composition is not encoded -i.e. when elements and the reference species of the 
linear combinations method coincide- results converge and correspond to the 
same Legendre transformation. 
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The advantage of resorting to Legendre transformations lies on the fact that, 
during the demonstration, the physical sense of the variables remains explicit. 
This allows for an interpretation of the Lagrange auxiliary function, which has a 
clear physical meaning: it is the complete Legendre transform of the Gibbs 
energy with respect to the element composition variables, and its value is zero 
for equilibrium states. 

The transformed composition variables are those linear combinations of such 
variables independent of the degrees of advance. They are also the conjugate 
variables of element potentials. These potentials are: (1) the intensive variables 
which go with certain mole number combinations (appendix A); (2) the Lagrange 
multipliers when minimizing the Gibbs energy subject to stoichiometric 
constraints; and (3) the conjugate of element composition variables when 
obtaining the complete Legendre transform of G. 
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APPENDIX A 

TRANSFORMED COMPOSITION VARIABLES AND DEGREES OF ADVANCE 

Taking as a basis the system approached in section 2, here it will be demonstrated that 
the transformed composition variables are the linear combinations of the composition 
variables that turn out to be degree of advance-independent. 

By separating from the initial quantities of B and C the ones participating in each 
reaction, we obtain 
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n 1 77 01 

n c -n c 
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n B n B 
1 


n} c = 


itf-ft; 






n A = 


nl + fr; 






n\ = 


< + ft; 






A, ~ 


< - nf 


«D-«D 


„2 02 

n c -n c 



-1 -1 



Trans. Luv. Acad. 2002, 1, 23-36 



34 



J. J. Luetich 



Tin = n D -^2', 
n% = n° c 2 + 2fr. 

And, choosing adequate linear combinations of the composition variables, 

n A + n c + 2n B = n A + £. + n^ -% x +n c 2 + 2& + 2 (n B - &); 

= n A + n c + 2n D ; 
n B + n c + n D = n B + £1 + nf - & + n c x - ft + nf + 2 ft + n D - ft; 
= nl + n° c + n B . 
No £, corresponds to these variables: they do not advance. 

Note that the key of this reasoning is the separation of the initial quantities of B and C 
into two parts whose extensions are not arbitrary. 

APPENDIX B 

COMPARISON OF TWO TRANSFORMATIONS 

Transformations of composition variables in the Lagrange multipliers and the Legendre 
transforms techniques is previous to the application of these methods and can be written 
as follows: 



n A \ 
n B 
n c 
n B ) 



n , 



n\ 



On the contrary, the transformation of the linear combinations technique is the result of 
the application of this technique and can be expressed as follows: 



n A 
n B 



B 



n B 



n A 
h B 



The circumstances in which they lead to the same result can be found solving the 
equation 
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The components of C are 
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Cll C12 C13 C14 

C21 C 22 C23 C 2 4 



1 -6n -6 12 

1 -621 -&22 



Thus, they lead to the same result when matrix C can be divided into a diagonal matrix 
and another one that is the opposite of B, i.e. when elements are selected among the set 
of components. 

APPENDIX C 

DIFFERENTIAL AND INTEGRAL EXPRESSIONS IN THE LEGENDRE 
TRANSFORMATIONS TECHNIQUE 

In section 2, the problem was solved by means of the Legendre transformations method, 
using differential equations. The same problem can be solved using these equations in 
their integral form. In this case, the method becomes the inverse of that introduced by 
Barbosa and Doherty (1987): given equations to obtain n* from n, the goal is to obtain 
the relation between £ and ju . 



G 



iU &) 



10 12 
111 



n B 
n c 

V ^D ) 



= 



(fiA-U /*B-& yUC-^A-^B /iD~2&-£ B ) 



n B 
n c 

V ^D ) 







This condition must be fulfilled when not all mole numbers are null (trivial solution), 

( yU A jU B yU C yU D ) = ( ^A ^B U + £} 2 £ A + ^ B )• 
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